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Abstract 

In this paper, we investigate the Cauchy problem for the compressible nematic liquid 
crystal flows in three-dimensional whole space. First of all, we establish the time decay 
rates for compressible nematic liquid crystal flows by the method of spectral analysis and 
energy estimates. Furthermore, we enhance the convergence rates for the higher-order 
spatial derivatives of density, velocity and director. Finally, the time decay rates of mixed 
space-time derivatives of solution are also established. 
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1 Introduction 

In this paper, we investigate the motion of compressible nematic liqnid crystal flows, which 
are governed by the following simplified version of the Ericksen-Leslie eqnations 

{ pt + div{pu) = 0, 

{pu)t + div(pu 0 -u) — pAu — {p + i^)VdivM -|- VP(p) = —yVd • Ad, (1.1) 

dt + u ■ Vd = 9{Ad + \Vd\‘^d), 

where p,u,P{p) and d stand for the density, velocity, pressnre and macroscopic average of the 
nematic liqnid crystal orientation field respectively. The constants p and z/ are shear viscosity 
and the bnlk viscosity coefficients of the fluid, respectively, that satisfy the physical assumptions 

/i > 0, 2p + ?)V >d. (1.2) 
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The positive constants 7 and 9 represent the competition between the kinetic energy and the 
potential energy, and the microscopic elastic relaxation time for the molecular orientation field, 
respectively. For the sake of simplicity, we set the constants 7 , 9 and ^*'(1) to be 1. The symbol 
0 denotes the Kronecker tensor product such that u®u= To complete the system 

fll.ll) . the initial data is given by 

ip,u,d){x,t)\^^Q = {po{x),Uo{x),do{x)). (1.3) 

As the space variable tends to inhnity, we assume 


lim {po - l,uo,do - wq)(x) =0, 

\x\—^oo 


(1.4) 


where Wq is an unit constant vector. The system is a coupling between the compressible Navier- 
Stokes equations and a transported heat flow of harmonic maps into S'^. It is a macroscopic 
continuum description of the evolution for the liquid crystals of nematic type under the influence 
of both the flow held u and the macroscopic description of the microscopic orientation conhgu- 
ration d of rod-like liquid crystals. Generally speaking, the system fll.ip can not be obtained any 
better results than the compressible Navier-Stokes equations. 

The hydrodynamic theory of liquid crystals in the nematic case has been established by Erick- 
sen and Leslie during the period of 1958 through 1968. Since then, the mathematical theory 
is still progressing and the study of the full Ericksen-Leslie model presents relevant mathematical 
difficulties. The pioneering work comes from Lin and his partners For example, Lin and 

Liu obtained the global weak and smooth solutions for the Ginzburg-Landau approximation 
to relax the nonlinear constraint d G They also discussed the uniqueness and some stability 
properties of the system. Later, the decay rates for this approximate system are given by Wu | 8 | 
for a bounded domain and Dai j^, 10 1 for the Gauchy problem respectively. Recently, Liu and 
Zhang 111], for the density dependent model, obtained the global weak solutions in dimension 
three with the initial density po ^ which was improved by Jiang and Tan 1^ for the case 


Po £ > §)• Under the constraint d G 5^, Wen and Ding [l3| established the local existence 


for the strong solution and obtained the global solution under the assumptions of small energy 
and positive initial density. Later, Hong 1^ and Lin, Lin and Wang ISjj showed independently 
the global existence of a weak solution in two-dimensional space. Recently, Wang [l^ established 
a global well-posedness theory for rough initial data provided that ||mo||bmo-i + [dolsMO < for 
some £0 > 0- Under this condition, Du and Wang [l7| obtained arbitrary space-time regularity 
for the Koch and Tataru t^e solution (u, d). As a corollary, they also got the decay rates. Very 
Recently, Lin and Wang [l 8 | established the global existence of a weak solution for the initial¬ 
boundary value or the Gauchy problem by restrict!^ the initial director held on the unit upper 
hemisphere. For more results, readers can refer to 19l-l24l| and references therein. 

Gonsidering the compressible nematic liquid crystal hows fll.ip . Ding, Lin, Wang and Wen 


25| have gained both existence and uniqueness of global strong solution for the one dimensional 


26 


space. And this result about the classical solution was improved by Ding, Wang and Wen 
by generalizing the huids to be of vacuum. For the case of multi-dimensional space, Jiang, 


Jiang and Wang |27| established the global existence of weak solutions to the initial-boundary 
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problem with large initial energy and without any smallness condition on the initial density 
and velocity if some component of initial direction field is small. Recently, Lin, Lai and Wang 


28l | established the existence of global weak solutions in three-dimensional space, provided the 


initial orientational director field do li^s in the hemisphere S^- Local existence of unique strong 
solution was proved if that the initial data (po, Wn, do) was sufficiently regular and satisfied a 


for the possible breakdown of such local strong solution at finite time could be found in 30 


The local existence and uniqueness of classical solution to (II.Ij) was established by Ma in |33|. 


30 

32 

in 

33 


On the other hand, Hu and Wu |50| obtained the existence and uniqueness of global strong 


solution in critical Besov spaces provided that the initial data was close to an equilibrium state 
(1,0, d) with a constant vector d G S"^. For more results, the readers can refer to 35| that have 
introduced some recent developments of analysis for hydrodynamic flow of nematic liquid crystal 
flows and references therein. 

If the director is an unit constant vector, then the compressible nematic liquid crystal flow 
fll.ip becomes the compressible Navier-Stokes equations. The convergence rates of solution for 
the compressible Navier-Stokes equations to the steady state has been investigated extensively 
since the first global existence of small solutions in (classical solutions) was improved by 


Matsumura and Nishida 3^. For the small initial perturbation belongs to only, Matsumura 


37( 1 took weighted energy method to show the optimal time decay rates 

\\V\p-l,u)it)\\L2<{l+t)--2 


for k = 1,2 and 


ll(p-1 ,m)(^)IU- < (1 + ^)' 


Furthermore, for the initial perturbation small in HH^, Matsumura and Nishida 38| obtained 

n with m > 4, Ponce 391 proved the 


L. <{l+t) 2^^ 2 


and for the small initial perturbation belongs to 
optimal decay rates 

for 2 < g < cx) and 0 < d < 2. With the help of the study of Green function, the optimal 
L'^(l < q < oo) decay rates were also obtained 40l-l42|| for the small initial perturbation to 
dd"* n with m > 4. These results were extended to the exterior problem (4^ or the half 


space problem [45|, 146| or with an external potential force |47j, but without the smallness of 
norm of the initial perturbation. While based on a differential inequality, Deckelnick in 


48l . l49j obtained a slower (than the optimal) decay rate for the problem in unbounded domains 
with external force through the pure energy method. Recently, Guo and Wang [l| developed a 
general energy method for proving the optimal time decay rates of the solutions in the whole 
space as 




2 < 
H^-l- ~ 


(1 t) 


-(l-i-s) 


(1,5) 


for 0 < / < iV — 1 by assuming the initial data ||(po — 1)'?^o)||ir-s(s G [0, |)) is finite additionally. 


This result was improved by Wang to establish the global existence of solution by assuming 


the smallness of initial data of dd^ rather than ddl- 2 . 


norm. 
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In this paper, we establish the global solution by the energy method [l(| under the assump¬ 
tion of smallness of ||(po — 1, ^o, V(io)||H3- The difficulty for the system (11.Ih is to deal with the 
nonlinear term Vd ■ Ad and the supercritical nonlinear term | Vdpd. Since this proof is standard, 
we only sketch it as an appendix for brevity. Furthermore, by assuming that ||(io ~ and 

II (po — l,Mo,do — tco)||Li are finite additionally, we establish the time decay rates for the com¬ 
pressible nematic liquid crystal flows by the method of spectral analysis and energy estimates. 
In order to improve the time decay rates of higher-order spatial derivatives of solution, we take 
the strategy of induction that is promised by the Fourier splitting method. Finally, we also study 
the decay rates for the mixed space-time derivatives of density, velocity and director. 

Notation: In this paper, we use G M) to denote the usual Sobolev spaces with 

norm || ■ ||ji/s and < p < oo) to denote the usual spaces with norm || ■ H^p. The 

symbol V* with an integer I > 0 stands for the usual any spatial derivatives of order 1. When I 
is not an integer, stands for defined by A^/ := f), where ^ is the usual Fourier 

transform operator and its inverse. We also denote <^{f) '■= /. The notation a means 
that a < Cb for a universal constant C > 0 independent of time t. For the sake of simplicity, we 
write ||(A,5)||x = \\A\\x + \\B\\x and / fdx := fdx. 

Now, we state our hrst result concerning the global existence of solution to the compressible 
nematic liquid crystal flows fll.lll - fll.Sp as follows. 

Theorem 1.1. Assume that (po — 1 ,mo, Vdo) £ for any integer N > 3 and there exists a 
constant <5 > 0 such that 

IKPo - l,Mo, Vdo)||H3 < d, (1-6) 

then the problem fll.ip - fll.3p admits a unique global solution (p, u, d) satisfying for all t >0, 

||(p- 1 ,W, Vd)||^. + A||Vp||^.-i + \\{Vu,V^d)\\lr.)dT < C\\{po - l,Wo, Vdo)||^.. (1.7) 

Jo 

After having the global existence of solution for the compressible nematic liquid crystal flows 
fll.ll) - fll.3p at hand, we hope to investigate the long-time behavior of solution. Hence, we establish 
the following time decay rates. 


Theorem 1.2. Under all the assumptions of Theorem M.li assuming the initial data ||do —tColU^ 
and IKpo — l,Mo,do — tco)||Li cire finite additionally, then the global solution {p,u,d) of problem 
fll.ll) - fll.4l) satisfies 

||V‘(p - l)(i)ll«»- + ||V‘«(i)llH»-» < C(1 +1)-^. 

where k = 0,1 ,..., iV — 1 and / = 0 , 1 , 2 ,..., iV -|- 1 . 


Remark 1.1. For any 2 < p < 6, by virtue of Theorem \1.2\ and Sobolev interpolation inequality, 
we also obtain the following time decay rates: 


V\p - l)(f)|Up + ||V'=M(f)|Up < ^(1 + 
V\d - wo)mL^ < C{1 + 
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where k = 0, — 2 and I = 0,1, 2,iV. Furthermore, it is easy to get the convergence 

rates, 

||v'=(p- i)(t)|Uoc + ||v'^M(t)|U.o < ^(1 +1)-^, 

||V*((i — u'o)(t)||L°° ^ ^*(1 Ft) ^ , 
where /c = 0,1, N — 3 and / = 0,1, 2..., N — 1. 


Remark 1.2. Generally speaking, the application of spectral analysis and energy estimates only 
helps the higher-order spatial derivatives of solution obtain the same decay rates as the first-order 
spatial derivatives of solution when we deal with the compressible Navier-Stokes eguations. In 
order to imwove the convergence rates for the higher-order spatial derivatives of solution, Guo 
and Wang Jj/ developed a general energy method when the initial data belongs to some negative 
Sobolev space additionally. In this paper, we apply the Fourier splitting method by Schonbek 15^/ 
or 153/ to improve the convergence rates for the higher-order spatial derivatives of solution. The 
advantage of our results fll.Sp is to verify that fll.Sp holds on for the case s = |. Furthermore, 
the time decay rates in m 


jN+l 


id-wo)mL^<cii+t)-^, 


is completely new. 

Finally, we turn to study the last results concerning the time decay rates for the mixed 
space-time derivatives to the compressible nematic liquid crystal flows fll.ip - fll.4p . 

Theorem 1.3. Under all the assumptions of Theorem \1.2i then the global solution {p,u,d) of 
problem ffTTD-ffLll) has the time decay rates: 

IIVViWIlr.™ + l|V^*(f)|U= < d/(l + t)-^, 

||V'di(f)|U2<C'(l + f)-^ 
for k = 0,1,..., N — 2 and / = 0,1,..., iV — 1. 

This paper is organized as follows. Since the proof of global existence of solution is standard, 
we only sketch it in section 4 as an appendix. In section 2, we establish the time convergence 
rates for the density, velocity and direction held to the compressible nematic liquid crystal hows 
fll.ll) - fll.4p . More precisely, the decay rates are built by the method of spectral analysis, energy 
estimates and Fourier splitting method. In section 3, we also study the time decay rates for 
mixed space-time derivatives of density, velocity and director. 

2 Proof of Theorem 

In this section, we investigate the time decay rates for the compressible nematic liquid crystal 
hows fll.ll) - fll.4p when the initial data belongs to space additionally. First of all, we derive 
the decay rates for the linearized compressible nematic liquid crystal hows that are a coupling of 
linearized compressible Navier-Stokes equations and heat equations. Secondly, we establish the 


1.2 
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decay rates for the compressible nematic liquid crystal flows fll.ip - fll.4l) by the method of spectral 
analysis and energy estimates. Furthermore, we improve the decay rates for higher-order spatial 
derivatives of density, velocity and direction field. 

Denoting g = p — 1 and n = d — Wq, then we rewrite fll.ip in the perturbation form as 


Qt + divM = Si , 

ut — pAu — {p + z/)VdivM -|- Vp = 5*2, 
rit — An = S 3 . 

Here Si{i = 1, 2, 3) are defined as 

51 = —gdivu — u - Vg, 

5 2 = —u ■ Vu — h{g)[pAu + {p + z/)VdivM] — f{g)'Vg — g{g)'Vn ■ An, 
S's = —M ■ Vn |Vn|^(n wq), 

where the three nonlinear functions of g are defined by 

Q 


( 2 , 1 ) 


( 2 , 2 ) 


h{Q) - 


fie) '= ^ - 1 and g(e) := 


(2.3) 


(2,4) 


P+1 P+1 P+1 

The associated initial condition is given by 

{g,u,n)\^^Q = {go,uo,no). 

By virtue of (11.711 and Sobolev inequality, it is easy to get 

2 - - 2 

Hence, we immediately have 

l^(p)U/(p)l < C'IpI and |p''"^(p)|, |h''(p)|, |/^(p)| < C for any k>l, (2.5) 

which we will use frequently to derive the a priori estimates for the time decay rates. 

Now, we state the classical Sobolev interpolation of the Gagliardo-Nirenberg inequality, refer 
to 


54 


Lemma 2.1. Let 0 < ni,a < I and the function f G then we have 

where 0 < 6 * < 1 and a satisfy 


( 2 , 6 ) 


a 


m 




I 


P 

On the other hand, the following lemma is very useful when we deal with the nonlinear 
function of p, refer to M • 

Lemma 2.2. Assume that ||p||loo < 1 . Let g{g) be a smooth function of g with bounded deriva¬ 


tives of any order, then for any integer m > 1 we have 


|V"^(p(p))IU^ < ||V"*p|| 


L°°- 


(2,7) 
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2.1. Decay rates for the nonlinear systems 

First of all, let us to consider the following linearized compressible nematic liquid crystal 
systems 

{ Qt + divM = 0, 

Ut — /iAn — (p + i^)Vdivn + Vp = 0, (2.8) 

Ut — An = 0, 

with the initial data 

ig,u,n)\^^f^ = {go,Uo,no). (2.9) 

Obviously, the solution {g,u,n) of the linear problem fl2.8l) - fl2.9p can be expressed as 


{g, u, n)*’’ = G{t) * {go, uq, no)*^ t > 0. (2.10) 

Here G{t) := G{x,t) is the Green matrix for the system (12.81) and the exact expression of the 
Fourier transform G{f,t) of Green function G{x,t) as 


G{f,t) = 


where 


Since the systems 


A+e^-*-A_e^+* 




A+e^+^-A-e^-‘ (j _ 

A+-A- A+-A_ [Dx3 |5|2^ 

0 


A+-A_ 

f 
0 


Ao = -/i|eP, Ai = -|ep, 






0 

0 

e"i*/3x3 


is a decoupled system of the classical linearized Navier-Stokes equations 


and heat equations, the representation of Green function G{f,t) is easy to verify. Furthermore, 


we have the following decay rates for the linearized systems fl2.8p - fl2.9l) . refer to 50 


Proposition 2.3. Let N > 3 be an integer. Assume that {g,u,n) is the solution of the linearized 
compressible nematic liquid crystal system fl2.8p - fl2.9p with the initial data {go,Uo,no) € H^DL^, 
then 

IIvVlIi. < C (||(e„, «„)|||. + IIv‘(e„,«„)||i.) (1 + 

IV'-Miii. < c + ||v'-(£>o,t.o)lli.) (1 + 

II V‘n||i, < C (llnolli. + llV^nolll,) (1 + «)-§-" 

for 0 < A; < A. 


The following estimates are essential for us to establish the time decay rates by the method 
of Green function. Since it is easy to derive, then we only state the results here for the brevity. 
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To be precise, we have 

ll(Si.%,S3)||i. < (llelli. + I|t.||i. + llVnlU^XIIVelU^ + ||Vt.||„. + ||Vn||„.) 

<i5(l|Ve|U« + ||Vti||„.+ ||Vn||„.), 

II('S'l, S'2, 53)11^2 < (llplli^i + llwlliyi + II Vn||jyi)(|| V^p|| 2,2 + ||V^m||hi + ||V^n||j^i) 

<5(||V2p|U2 + ||V2M||Hi + ||V2n||^,i), 

||V(^i,^2,^3)||l 2 < (||p||^,2 + ||m||^, 2 + ||Vn||^,2)(||V2p||^^i + IIV^mIIhi + ||Vnb2) 
<5(||V2pbi + ||V2M||^,i + ||Vn||^,2). 

Before studying the time decay rates, the following energy estimates will be used to guarantee 
the hrst-order derivatives of velocity and director enjoying the same convergence rates. 


Lemma 2.4. Under the assumption fll.61) . then we have for any integer k = 0,1,2,N, 

d 


dt 


\V^n\^dx+ / iV^+^nl^da; < (5||V^+^n||i2. 


( 2 . 12 ) 


Proof. Taking k—th spatial derivatives to fl2.ip n. multiplying the resulting equations by and 
integrating over M^(by part), it arrives at 






= / [-U ■ Vn + |Vnp(n + wq)] • V^n dx = h + h. 


(2.13) 


Since the inequality (I2.12p is easy to verify for the case k = 0, then we only verify the case k>l. 
Applying Leibnitz formula and Holder inequality, it is easy to deduce 

k-l 


h = 


/ ET-i 

1=0 


VuV^-^n dx 


Al + l, 


k-l 


(2.14) 


<5^l|V'«|U,||V‘-'n||i.||V‘+‘n||i3. 


1=0 


For the case 0 < / < , by virtue of the interpolation inequality fl2.6p and Young inequality, 

we have 

II Y7^+l/> 

L6 Y 


|V^n||r3||V^"'n| 


< II 


u\ 


Il 2 I 


V'=+'n||| 2 ||Vn||| 2 ||V''+V||;^ 2 "||V''+'n|U 2 


n\\L^ 

I 


fc+i„ 


k |IV7^+l,. 


(2.15) 


< 

rsj 


where a is dehned 


<5(||V^+'n||i2 + ||V^+V||i 2 ), 

a=l 


0 ,- 

’ 2 


2{k-l) 

Similarly, for the case + l</</c — 1, it follows that 

i 6 ||V'^+V|L 2 


|V'n|U3||V^“'n| 


1 +^ 




1 +^ 


^ IIWIII2IIV^+^n|| 2 r IIV“n|| 2 r II|| 


Ll 2 


< 

rsj 


5(||V'^+in||i2 + ||V'=+^n||^ 2 ), 




(2.16) 
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where a is defined by 


a = 


k + l 
2 / + 1 




Indeed, the smallness of ||V“n||(0 < a < 1) is guaranteed by the boundedness of ||ri||L 2 and the 
smallness of ||Vn||i 2 . Plugging fl2.15p and fl2.16p into fl2.14p . we obtain 

/i<h(||V'=+'M||i. + ||V'=+in||i.). 

By virtue of Leibnitz formula and Holder inequality, it arrives at 


(2.17) 


h = - 


V^-WVn\'\n + m;o))V^+V dx 

k-l 


/=0 

/ k—1 I 

EE + wo)V’^^^n dx 

1=0 m=0 

/ » fc— 1 

I VnpV^-^nV^+^n dx - C'r-iV”^+^nV^-”*n(n + wo)V’^+^n dx 

'' m=0 

k-2 I 

-/EE 

Z=1 m=0 


(2,18) 


(2.19) 


n dx 


— hi + h2 + hs- 

With the help of Holder inequality and interpolation inequality fl2.6p . we deduce directly 
hi < ||Vn|U6||Vn|U6||V'^-in|U6||V'^+in|U2 

< ||V2ri|U2||Vn||^;^||V"+'ri|||,||Vn|||2||V"+'ri||i^;^||V"+V|U2 

< ||V 2 n|U 2 ||Vn|U 2 ||V"+'n||i. 

In order to estimate the term h 2 , using Holder inequality and fl2.6l) . we obtain, for the case 
0 < m < [^], that 

||V"*+^n||i3||V^-™n||i6||V^+^n||i2 

-f_m m m -i_m 

< II V“n||y ‘ ||V‘+‘n||^‘,||Vn||^‘,||V‘+‘n|i;,; ‘ II V‘+'n||i2 

< illV^+'nlli., 

where a is defined by 


( 2 . 20 ) 


CK — 1 -|- 


k 




2{k — m) 

Similarly, for the case -|-l<m</c — 1, itis easy to deduce 

||V"^+^n|L3||V^-™n|L6||V^+^n|L2 


rn-f--^ rto-r-^ 

< ||Vn||^~||V'=+'n||;;;;^||V“n| 


< 

n>j 


lL 2 

5||V^+'n||i„ 


1,2 


2 

fc iiy 


1 "*■'■2 
^+i,>nlr fc livyfc+i 


n\ 


L 2 


n\\L^ 


( 2 . 21 ) 
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where a is defined by 


O' — 1 -|- 


k 


2 m + 1 


Combining fl2.20p with fl2.2ip . then the term I 22 can be estimated as follow 

( 2 . 22 ) 

To deal with the term J23. For the case 1 < / < by fl2.6p and Holder inequality, we have 


V”^+^n|b6||V 


/+!—m 


1 -- 




l — Z —m+1 

< Iiv“nii;3^iiv"+vii72iivnii;;^iiv"+'^ 


lL2 
Z + 1 


n\ 


L2 


X ||Vn||,,S l|V‘+‘n||l; 


, 1 - 


1+1 




1 m 1 I m 

k llY7.>nll + k ||Y7fc+l.„||2 


< II V"n ||^2 II Vn|| 2 I ^ IIV'-' "n 


V '^n\\L 2 

L2 


< 




where a is defined by 


(y. — 1 T 


k 


k — m 


e|2,3). 


Similarly, for the case + 1 < / < fc — 2, it is easy to obtain 

||V"*+^n|b6||V'+^“™n|b6||V^"^“'n||i6||V^+^n|b2 


< 


liVnII 


1- 


m+1 


jk+l. 


m +1 


liVnII 

1+2 , ^ ^ 1+2 


1,2 II V n \\^2 


-| I —m+1 I —m+1 

k ||V7fc+l.nll f'' 

L2 II V ^11^2 


X V“n 


“.nil k llvyfc+lnlr k 


L2 


n\ 


k+1 


L2 


n\\L^ 


<||Vn||^;'^'||V“n||J||V^+' 


< 


lL2 

<5||V^+^n||i2, 


n\\l 2 


where a is defined by 


a — 1 T 


k 


1 + 2 


€ [2,3). 


Combining fl2.23p and fl2.24p . it follows that 

/23 < ^llV^+'nlli., 

Substituting fl2.19p . fl2.22p and (12.251) into fl2.18p . then we get 

/2<3||V‘+‘n|||2. 

which, together with fl2.17p . completes the proof of lemma. 


(2.23) 


(2.24) 


(2.25) 

(2.26) 
□ 


Now, we turn to establish the time decay rates for the compressible nematic liquid crystal 
flows fi2.ip - fi2.4p . 
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Lemma 2.5. Under the assumptions of Theorem \1.2i the global solution {g,u,n) of problem 
(12.II) - (12.41) satisfies 

II + ||V"n(t)||^^+i_. < C(1 + 

for k = 0 , 1 . 

Proof. First of all, choosing the integer k = I & [0; -^] Lemma [2.41 we have 


(2.27) 


fL 

dt 


V'nli:, + C llV'+^nli:, < 6 ||V'+^m 


Il 2 


L 2 ~ 


The combination of fl4.7p and fl2.28p yields 
d 


+ Cr (II+ IIV'+^w||l,„ 



(2.28) 

+ II V^~*~^?7.||^m+l-i) Y 0, 

(2.29) 


where iPpit) is dehned as 

jr-(t) = ||v'(^,m)||^_* + ||V'n||^^+i_, + 


i„l|2 , 2C'2(5 y-v 


a. 


V'^u-V’^^^gdx. 


l<k<m—l ' 


With the help of Young inequality, it is easy to deduce 


a- 


-1 ( II V7i 


v'(p,m) 


I 


IV'n"^ 


fjm + l 




i^l|2 


(2.30) 


Adding on both sides of fl2.29p by ||V^(p, M,n)||^2 and applying the equivalent relation fl2.30p . 
then we have 

fTr{t) + GTr(t) < \\V\s,u,n)\\l,. (2.31) 

Taking I = 1 and m = N specially in (I2.3ip . it arrives at 


fL 

dt 


J^^{t) + CJ^f{t)<\\W{Q,u,n)\\i„ 


which, together with the Gronwall inequality, gives 


(t) < ^(^(0)e-^‘+ / e-^(*-")||V(p,w,n)||i,dr. 


(2.32) 


In order to derive the time decay rate for {t), we need to control the term ||V(p, m, n)||^2- In 
fact, by Duhamel principle, one represents the solution for the system (I2.ip - (l2.4p as 


{Q,u,nY^{t) = G{t)*{Qo,uo,nor + / G{t - s) * {S^, S 2 , {s)ds. 


(2.33) 


Denoting F{t) = sup (1 + r)2(|| Vp(r)||^jv-i + II VM(r)||^^_i -f ||Vn(r)||^jv), by virtue of (|2 .11 | ) , 

0<T<t 

fl2.33p and Proposition 12.31 then we have 

II V(p, w, n)||i2 < G(1 + f)-i + G r (II (Fi, ^2, ^3)|lii + || V(Fi, F2, F3)||i2) (1 + t - r)-idr 

Jo 

< G(1 + t)-l +G f 6 (II Vp||^2 + II Vw||^2 + II Vn||^2) (1 + f - r)-tdr 


<G(l + f) 2 + G( 5 F(f) / ( 1 +t —r) 2(1 + 7-) 2dr 

Jo 

< G{l + t)-l +G 6 F{t){l + t)-L 
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where we have used the fact 

rt 


(1 + t — r) '^{1 + t) ^dr 


10 


+ / ( 1 +t —r) 2(1 + 7-) ^dr 


'0 


< {i + tyK 

Thus, we have the estimate 

\\Vig,u,n)\\y<Cil + t)-l{l + 6Fit)). 

Inserting (12.34^ into (12.321) . it follows immediately 

< F^{0)e-^* + C [ + r)"t(l + 5F{T))dT 

Jo 

<F^{0)e-^^ + C{l + 6F{t)) f e-^(*-^)(l-f r)-idr 

Jo 

< F^{0)e-^^ + ^(1 + SF{t)){l + t)-l 

< C{l + 6F{t)){l + t)-l, 

where we have used the simple fact 

[ e-^('-")(l + r)"5dr 

Jo 

= [ + [ + T)~^dT 


< e 2* / (1 + r) '^dr + [1 + - 




< C(l + t)"2 . 

Hence, by virtue of the definition of F{t) and (12.351) . we have 

Fit) < Cil + 5Fit)), 

which, in view of the smallness of <5, gives 

Fit) < a 

Therefore, we have the following decay rates 

||Vp||^jv-i + llVull^AT-i + ||Vn||^Ar < (7(1 Ft) 2. 


(2.34) 


(2.35) 


(2.36) 
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On the other hand, by fl2.1ip . fl2.33p . fl2.36p and Proposition 12.31 it is easy to deduce 

< 0(1 + t) 2-1-0/ ^11(^5'^^ 5*2^ 5*3)||2^-)-11(^5'^^ 5*2^ 5*3)1122^ (^1-l-1 _-r) 2 (ir 

Jo 

< 0(1 + t) 2-|-O / ^ (II Vpll^i + II Vm ||^2 + II Vnjl^s) (1 + t — t) ^dr 

Jo 

<C{l+t)-l + C [ (l + t-r)-i(l + r)-idr 

Jo 

< 0 (l + t)-i, 

where we have used the fact 

(1 + t - r)“5(i + r)“^dr < O (1 + t)~^ . 

Hence, we obtain the following time decay rates 

||(p,«,n)||i2 <o(i + f)-i. 

Therefore, we complete the proof of the lemma. □ 

2.2. Improvement decay rates for the higher-order spatial derivatives of solution 

In this subsection, one will improve the time decay rates for the higher-order spatial derivatives 
of density, velocity and director. More precisely, we have the following convergence rates. 

Lemma 2.6. Under the assumptions of Theorem \1.2\. the global solution {g,u,n) of problem 
(I2.1jl - (l2.4j) has following decay rates for all t > to(to is a constant defined below), 

II V^^ll^iv-. + + IIV'^nll^^+i-, < 0(1 + t)-!-’^ (2.37) 

where k = 0,1, 2,..., iV — 1. 

Proof. We will take the strategy of induction to give the proof for the convergence rates (I2.37p . 
In fact, the inequality fl2.27p implies fl2.37p for the case k = 1. By the general step of induction, 

assume that the decay rates (12.371) hold on for the case k = I, i.e. 

||V^p||^Ar-i + II V^MlI^iv-i + II V*Ti||^iv+l-! ^ 0(1 -ft) 2 f (2.38) 

for / = 1,2,3, ...,N — 2. Then, we need to verify that (I2.37P holds on for the case k = I 1. 

Indeed, replacing / as / -|- 1 and taking m = iV in fl2.29p . we have 

^J-z^i(t) + O 7 (IIV'+2^||^._._2 + II < 0, 

which implies 

+ Y + l|V'+'^?||^iv-*-2 + + ||V'+"n||^^_,) < 0. (2.39) 
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Denoting the time sphere S'o(see |^) as follows 


5'o:=UeM'lie|< 


R 


1 + t 


where i? is a constant defined below. By virtue of Parseval identity, then it is easy to deduce 

Jr3/So 

R f 


> 


> 


^ + t Ju3/Sq 

R 




> 


1 + t 
R 

1 +1 






(l + «p 


So 


leri^ix 


R^ 




a + ty 


Hence, we have the following inequality 




I „l|2 


1 + t 


Similarly, it is easy to obtain 




R 




a + ty 


R^ 


IV'^II 


L 2 - 


1 + t"' {i + ty' 

Summing up in fl2.4ip with respect to k from k = l to k = N — one deduces 






lL 2 - 






iH— - YT~t 
In the same manner, it is easy to deduce 


IJN-I-I 


{^+ty 




ffN-l-l. 


l+l^\\2 




R^ 


\ hn-i 


(1 + t) 


2 


Substituting (I2.40p . (12.421) and (I2.43p into (I2.39p . it follows immediately 

R 




^ (||V'+V|li2 + + ||V'+'n||^._,) 


< 


CjR^ 


1 +1 

V*p||^2 + II V^wll^jv-l-i + ll^^^ll^JV-i) 


2 (l + t )2 

For some sufficiently large time t > i? — 1, we have 

R 


1 + f 


< 1 , 


(2.40) 


(2.41) 


(2.42) 


(2.43) 


(2.44) 
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which implies 


R 




1+2 \\2 


1 + t 


(2.45) 


Plugging fl2.45l) into (12.441) . it arrives at 


J + 2(1 + t) 


< 


R‘^C7 


+ IIV^nll^jv-i-i + IIV^nll^jv-i) , 


2(l + t)2 

which, together with the equivalent relation fl2.30p and the convergence rates fl2.38p . yields 




-t;-l 


(2.46) 


Choosing 


R = 


2(/ + 3)C*g 

cv ’ 


and multiplying both sides of fl2.46p by (1 + we have 

I [(l + O'^T-,!,(()] <C(l + ()-i, 

for any t>to and to := — 1. Integrating fl2.47p over [0,f], it follows directly 

< h+i(0) + C(1 + t)^ (1 + 


(2.47) 


which, together with equivalent relation fl2.30p . gives rise to 

IIV'+Vll^iv-*-! + IIV'+^n||^^-*-i + IIV'+V||^2v-* < 0(1 + f)-i-'. 

Hence, we have verified that fl2.37p holds on for the case k = / + 1. By the general step of 
induction, we complete the proof of the lemma. □ 

Now, we will focus on improving the time decay rates for the iV—th and {N + 1)—th order 
derivatives of direction field. More precisely, we will establish the following time decay rates. 

Lemma 2.7. Under the assumptions of Theorem, \1.2i the global solution {g,u,n) of problem 
(12.1|) - (12.41) satisfies 

II < ^(1 + t)-!-’^ (2.48) 

for integer k = 0,1, 2 ,..., iV -|- 1 . 

Proof. Taking TV— th spatial derivatives to both sides of ( 2 . 1 ) 3 , multiplying V^n and integrating 
over M^, then we have 


1 d 

2 dt 







j V^(m ■ Vn)V^ndx + 


j V^(|Vn|2(n + woW^ndx. 


(2.49) 
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By virtue of the integration by part, Leibnitz formula, Holder and Young inequalities, it is easy 
to deduce 

- j V^(n • Vn)V^n dx 

Af-l 

N-I (2.50) 

^ \ ^ IIV7^,i.ll II II II V7-^+l/>n II _L lU.II II II 2 

~ / II Y V ?T.||£2||V ^||l2 + ||M||j|/i|| V ^||_L2 

k=l 

N-1 

^ E l|V‘«irff.||V~'"'"‘«lli- + (<r + <5)||V‘^+‘r!||i.. 

k=l 

Taking k = N in fl2.18p specially, we have 

j V^(|Vn|2(n + n;o))V^nda; < S\\V^+^n\\l 2 . (2.51) 

Substituting fl2.50p and fl2.5ip into fl2.49l) . in view of the smallness of 6 and e, it arrives at 

N-l 

^dx < II (2.52) 

k=l 


d 

dt 




|V^+V| 


On the other hand, taking {N + 1)—th spatial derivatives to both sides of (2.1)3, multiplying by 
and integrating over M^, then we have 


lA f I\V^^^n\^dx 


= -J V^+^(n ■ Vn)V^+^nda; + J V^+^(|Vn|2(n + M;o))V^+^n dx. 

By virtue of the Leibnitz formula. Holder and Sobolev inequalities, it is easy to deduce 


(2.53) 


- J Vn)V'''+'n dx 

N-l 


^ ll^^“IUH|V^+^"^n|U6||V^+2n||i2 + ||Vn|Uoc||V^n|U2||V^+2n|U2 


k=0 

N-l 


^ IIV7^ II IIV7.^“I“2—^ II IIV7.^“t"2 II I II II IIV7.^~t"2 

^ / ||V 'u||/fi||V ?7.||i2 ||V ^||l 2 + ||'*^||hi|| V n 


N+2 ||2 

L 2 


(2.54) 


k=l 


< 

r\j 


+ ||Vn||Loc||V'"n|U2||V'"+"n||i2 
N-l 

k \\2 II y 7A^H-2—/c II2 


Y1 llV'^“ll^i||V'^^'"''^lli2 + ||Vn||ioo||V"n||i2 + (e + 6)\\V^- 


7 Af+ 2.„||2 


k=l 


On the other hand, taking /c = Y + 1 in (I2.18p specially, it is easy to deduce 

[ V^+^(|Vn|2(n + n;o))V^+V dx < <5|| 


(2.55) 
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Substituting fl2.54p and fl2.55p into fl2.53l) . by virtue of the smallness of 6 and e, we have 

d 






dt 


N-l 


< 


Y, IIV‘«lll,.||V"+''-'=n||i, + ||Vn|||„||V««||i,. 


(2.56) 


k=l 


Adding fl2.52p to fl2.56l] and applying the Sobolev interpolation inequality fl2.6l] . it arrives at 


^ [ (iV^np + iV^+^nHcix + f (|+ \V^+M^)dx 


N-l 


N-l 


< 


Y IV‘'“ll?,.||V"+‘-*n|ry + Y IIV‘«ll?,.||V''+=-‘n||=. 


k=l 


k=l 


+ \\Vn\\U\V^u\\l, 

= III + II2 + Ih. 

Applying the decay rates 02.371) . it follows directly 


v-i 


II, = l|V‘«||?,.||V''+‘-‘n||i. + ||Vu||'J,.||V"n|||, 


k=2 

N-l 


< ^(1 + + (1 + t)"i(l + 


k=2 


< (1 + + (1 + 


Similarly, we obtain 


N-l 


Ih < Y l|V‘«llt.||V"+^-‘n||i. + ||V«||y llV^+'nlli, 


k=2 

N-l 


< ^(1 + + (1 + 


k=2 


< (1 + + (1 + 
< (l+t)-(^+3), 


and 


Ih < (1 + < (1 + t)-(^+3). 

Inserting 02.58p - 02.60p into fl2.57p . it arrives at obviously 

^ j{\V^n\^+ \V^^^n\^)dx+ j+ \V^^^n\^)dx < {I + 


(2.57) 


(2.58) 


(2.59) 


(2.60) 


(2.61) 
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Taking the Fourier splitting method as the inequality fl2.40p and the decay rates fl2.37p . we have 

^ + + [{\V^n\^+\V^+\\^)dx 

at J l + t J 

~ + + (2.62) 

< (1 + + (1 + 

Multiplying fl2.62p by (1 and integrating the resulting inequality over [0, t], then it follows 

directly 

||V^n||^i < C(1 +t)“(i+^), 

which, together with (I2.37p . implies 

II < ^(1 + (2.63) 


for k = 0,1, 2,iV. On the other hand, from the inequality fl2.56p . it arrives at 

^ j\V^+\\'^dx+ j iV^+^nl^dx 

N-l 

< \\\7k,,\\2 


E llV'wll^illV^+'-'nlli^ + ||Vn||^i||V--^n 

k=2 
Af -1 

< ^(1 + ^yil+N-k) + (1 + t)-i(l 

k=2 

(l+^)-(V+5)^(l^^)-(V+4) 

(l^^)_(JV+4)_ 

Taking the Fourier splitting method as 02.401) . we have 

d 


Ih + ||VM||^i||V^+V|li2 + ||V^n||i2||V2n||^i 

+ t)-(|+^) + (1 + t)-(^+4) 


< 

rsj 

< 



^|V^+V|' 

1 + ^ J 



Htt? j 

< (1 + t)-(3+^) + (1 + t)-(^+T 


(2.64) 


(2.65) 


Multiplying 02.65P by (1 +f)'^’''^ and integrating the resulting inequality over [0, t], then it follows 
directly 

llV^+^nllis < C'(l + t)"(i+^), 

which, together with (I2.63p . completes the proof of lemma. □ 

Proof for Theorem 11.21 With the help of Lemma 12.61 and Lemma 12.71 we complete the 
proof of Theorem 11.21 
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3 Proof of Theorem 11.3 


In this section, we will establish the time decay rates for the mixed space-time derivatives of 
density, velocity and direction held. More precisely, we have the following decay rates. 


Lemma 3.1. Under the assumptions of Theorem \1.2[ the global solution {g,u,n) of problem 
fl2.lll - fl2.4ll has the time decay rates 


I < (1 + t) 2 


(3.1) 


for any integer k = 0,1,N — 2. 


Proof. Taking k—th (k=0,l,...,N-l) spatial derivatives to (12.ip ^. multiplying by and inte¬ 
grating over then we have 


II = - j V\dwu + gdiYu + u ■ VQ)V^Qtdx = Ilh + Hh + Hh- 

By virtue of the Young inequality, it follows directly 

///i<||v'=+‘«||i.+£||v‘a||i.. 

Applying the Leibnitz formula. Young and Sobolev inequalities, we have 

k 

iih < 5^I|v'^iU3||v"+'-'m|U6||v'=p 4|U2 + bl|Loc||v"+'M|U2||v"ft|U2 


1=1 

1=1 

k 




ioo||v''+^w|li 2 + £||vVdli= 

+ (l + t)-3(l+t)-i-'^ + £||V'=Pi 


k„ l|2 
L 2 


<5^(1+«)-"'(!+ i)- 

1=1 

<(i + t)-^-^ + £||v^p,lli- 

In the same manner, it follows immediately 

JJJ3<(l+t)-^-^ + £|lV^ft||i2. 

Inserting fl3.3p - fl3.5p into fl3.2p . then it is easy to deduce 

l|vVJi2<(i + t)-"-" + l|v'=+iM||i2 

which, together with the time decay rates fll.811 . completes the proof of lemma. 

Next, we establish the time decay rates for the space-time derivatives of velocity. 


(3.2) 


(3.3) 


(3.4) 


(3.5) 


□ 


Lemma 3.2. Under the assumptions of Theorem If.ill the global solution {g,u,n) of problem 
fl2.111 - 02.411 has the time decay rates 


|vV(()lli. <(i + ()-*-‘. 


(3,6) 


for any integer k = 0,1,..., N — 2. 


19 





















J.C.Gao, Q.Tao, Z.A.Yao 


Proof. Taking k—th (k=0,l,...,N-2) spatial derivatives to fl2.ip ^. multiplying by and inte¬ 
grating over then we have 


l|v‘u,|li.=/v‘(-u 


• Vu -|- g{g)[fiAu + {g. + i')'\/divu\)'\/’^utdx 


j V^[(/(^) + l)VQ+g{Q)VnAn]V^utdx 


= ivi + 1V2 + m + m- 


(3.7) 


Applying the Leibnitz formula, Holder and Sobolev inequalities, we have 


IV, < Y. IIV'«IU.||v‘+>-‘«iu.|| Will. 

1=0 

k 

1=0 

k 

< E(i + + <;|lvVlIi. 

1=0 

Similarly, we have 

k 

IV 2 < ||V'((7(^) - l)||Loc||V^+2-'M|U2||V^i|U2 + 

/=0 

k 

< E llv'elli-llv'-Wii^. + llvWili, + £||v%fy 

1=0 

k 

< E (1 + tp-‘0 + + liv‘+"«|li. + £||vVlli. 

1=0 

< (1 + + II vWlii, + ill Will.. 


(3.8) 


(3.9) 


and 


m < E ii^i/(rf + iiiu-iiv^+'-viii.iivviii. 

/=0 

k 

< EII v'ell!»llv‘+'-'e|l!. + Iiv^+Vlll. + £||v‘«,|||, 

/=0 

k 

< E(1 + «)■*■'(! + *)-«-“ + IIV^+Vlli. + ellvViii. 

1=0 

<(! + ()-"-'=+II v'=+V|li.+£||vVlli.. 


(3.10) 
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In the same manner, it arrives at 

/ k I 

EE 

1=0 m=0 
k I 

< ||V'^|Uoo||V”^+'n||^,i||V"+='-'-”^n|U2||V'^njL2 

Z=1 m=0 

+ \\g{g)\\L^\\Vn\\H4^>^^M\L4^%\\L^ 

1=1 m=0 

+ ll9(e)ll?,-l|Vn||J,.||V«n|| = . +£||V‘«,||i. 

k I 

< ^ ^(1 + t)-t-'(l + t)-t-”^(l + 

1=1 m=0 

+ (1 + t)-i(l + + e\\V’^u4l2 

<(l + t)-7-A3 + ^||Vfc«j2^. 

Inserting fl3.8p - fl3.lip into fl3.7l) and choosing e small enough, then we have 

\\^%it)\\h < (1 + + l|V'^+'M||i2 + II ^=+^1112. 

Therefore, we complete the proof of lemma. □ 

Finally, we establish the time decay rates for the space-time derivatives of director. 

Lemma 3.3. Under the assumptions of Theorem \1.2i the global solution {g,u,n) of problem 
fl2.II) - (12.41) has the time decay rates 

||V^n,(f)||i. <(l + ^)-i-^ (3.12) 


for any integer /c = 0,1,..., iV — 1. 

Proof. Taking k—i\i (k=0,l,...,N-l) spatial derivatives to fl2.ip n. multiplying by and inte¬ 
grating over M^, it arrives at 

V^(An -u-Vn+ \Vn\^{n + M;o))V^ntda; = l/i + I /2 + I/ 3 . (3.13) 



By virtue of the Young inequality, it follows immediately 

< IIV"+'n||i 2 + ellV"+'nt||i 2 < (1 + t)-l-^ + £||V"+'ni||i 2 . (3.14) 

Applying the Leibnitz formula. Holder and Young inequalities, we have 

^2 < 5^ l|V'n||^i||V'^+2-'n||i2 + e||V"ni||i2 
1=0 




-l-k+l 


+ e|| V^ni ||^2 


1=0 


< (l + f)-^-^ + £||V'=nt| 


1,2 • 


(3.15) 
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and 


fc-i i 


l|V”^+'ri|U6||V'+'-™n|U6||V"-'n|U6||V'^ni|U2 

1=0 m=0 
k 

+ ||V”^+'n||L3||V'+^-”^n|U6||V'=ni|U2 

m=0 

1=0 m=0 
k 

m =0 
1 / 

< ^ ^(1 + t)-i-'”(l + t)-i-^+^{l + t)-t-^+' 


m=0 
fc -1 I 


1=0 m=0 
k 


+ ^(1 + t)-t-™(l + + ellV'^nilli^ 


m=0 

^ /^1 I -/-^— 6 —/c 


(l + t)-^-'^ + £||V^ni "2 


L 2 - 


Substituting fl3.14p - fl3.16p into fl3.13p . it arrives at 

l|vMt)lli 2 <(i + t)-^^ 

Therefore, we complete the proof of the lemma. 


(3.16) 


□ 


Proof for Theorem 11.31 With the help of the Lemma 13.11 Lemma 13.21 and Lemma 13.31 we 
complete the proof of Theorem 11.31 


4 Appendix: Proof of Theorem 11.11 


In this section, we will establish the global existence of solution for the compressible nematic 
liquid crystal flows fll.ip - fll.3p . We only sketch it here since this proof is standard. More precisely, 
assume there exists a constant 5 > 0 such that 

\J ^o(^) = Il^(^)llrf3 + ll“(^)llr/3 + IIVn(t)||/^3 < 5. (4.1) 

Here, we denote g = p — 1 and n = d — Wo- Then, we derive some energy estimates that play 
an important role for establishing the global existence of solution under the assumption of fl4.1 1 . 
In fact, the following lemmas are easy to establish just following the idea by Guo and Wang [it]. 
Hence, we only state the results here for the sake of brevity. 


Lemma 4.1. If < 6, then we have 

^ |lV^(p, u, Vn)||", + G II Vn)||', < <51| ; 

I ||V'=+^(p, u, Vn)||", + G II Vn)||^, < 5 || ; 

^ y ■ V'^+^gdx + G|| V^+V|li2 < || V^+iM||i2 + || + || 
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where k = 0,1, 2,— 1. 


Proof for Theorem ll.lt Let N > 3 and 0 < / < m — 1 with 3 < m < N. Summing up the 
estimates (I4.2p ^ from fc = / to m — 1 , it is easy to deduce 


+Ci||V'+'(M,Vn)||^_*_i ^Cs^llV'+VlI^^-^-i- 

Taking k = m — 1 in (I4.2p „. then we get 

^||V”^(p,w,Vn)||i, + Ci||V"^+iKVn)||i, <C2<5||V>||i,. 

Adding 04.31) to 04.41) . it is easy to obtain 
d 


dt 




(4.3) 


(4.4) 


(4.5) 


Summing up the estimate 04.2l) o from k = I to m — 1, we obtain 

+ Csll V'+V||^^-<-i < C'4(||V' 

l<k<m—l 

Multiplying 04.6p by 2 C 25 /C 3 and adding to 04.51) . then we have 


+ ||V'+2Vn||^_,_i). (4.6) 


j^srifj+c, (||v'+‘^,|||,„_,.. + iiv‘+'«iii,, 


where Sp{t) is dehned as 


Sr{t) = \\V{g,u,Vn)\\l^.,+ 


2 C 26 




1-1 + — 0 ) 




l<k<m—l 


By virtue of the Young inequality, it is easy to obtain 

f I ■ V"+V| dx < 2 + \\Vu\\lrr.-t-i) . 

Hence, in view of the smallness of 6, we deduce 

Integrating 04.7p over [0,t], then we obtain 


(4.7) 


(4.8) 


£,"(*)+ft / (||V‘+Vlll,»-i-. + l|V‘+‘(t<,Vn)||^„_,)ciT<f,"'(0), 


l+ll 


which, together with 04.Sp . gives 


||V'(p(t),«(t), Vn(t))||^_, + / (||V'+V||^^-.-i + \\V^\u,Vn)rH^.,)dT 

Jo 

< C\\V^{go,Uo,Vno)\\Hm-i. 


(4.9) 
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Choosing I = 0 and m = 3 in fl4.9p . then we obtain 

MtWm + + ||Vn(t)||^3 < C{\\gorH^ + Wnorm + llVnoH^s). (4.10) 

By a standard continuity argument, the inequality fl4.10p will close the a priori estimate fl4.ip . 
Hence, taking I = 0 and m = iV in fl4.9p . it is easy to obtain 

||(p(t),M(t), Vn(t))||^iv + [ + II V(m, Vn) II j^]v)dT < CH (po) "Wo; Vno) || j^iv, 

Jo 

which completes the proof of Theorem 11.11 
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